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ABSTRACT 

Aims. We provided accurate estimates of distances, radii and iron abundances for four metal-rich Cepheids, namely V340 Ara, UZ Set, 
AV Sgr and VY Sgr. The main aim of this investigation is to constrain their pulsation properties and their location across the Galactic 
inner disk. 

Methods. We adopted new accurate NIR (J,H,K) light curves and new radial velocity measurements for the target Cepheids to 
determinate their distances and radii using the Baade-Wesselink technique. In particular, we adopted the most recent calibration of the 
IR surface brightness relation and of the projection factor. Moreover, we also provided accurate measurements of the iron abundance 
of the target Cepheids. 

Results. Current distance estimates agree within one cr with similar distances based either on empirical or on theoretical NIR Period- 
Luminosity relations. However, the uncertainties of the Baade-Wesselink distances are on average a factor of 3-4 smaller when 
compared with errors affecting other distance determinations. Mean Baade-Wesselink radii also agree at one cr level with Cepheid 
radii based either on empirical or on theoretical Period-Radius relations. Iron abundances are, within one cr, similar to the iron contents 
provided by Andrievsky and collaborators, thus confirming the super metal-rich nature of the target Cepheids. We also found that the 
luminosity amplitudes of classical Cepheids, at odds with RR Lyrae stars, do not show a clear correlation with the metal-content. This 
circumstantial evidence appears to be the consequence of the Hertzsprung progression together with the dependence of the topology 
of the instability strip on metallicity, evolutionary effects and binaries. 
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1. Introduction relations still lacks a firm empirical validation j[Benedict et al.l 

^ 2007; Romaniello et al. 2008; Sandage et al. 20^. 
. . ■ Classical Cepheids are used both as standard candle s and 

> ; tracers of young s tellar populations (Maciel & Costa. l2009t ™s thorny problem remains even if the zero-point and the 

■ lPedicellietal.ll2009l) . They are bright and variable objects and slope of the PL relation can be estimated with a variety of in- 

?S ' thanks to the Hubble Space Telescope they have been identi- dependent methods. Ideally the cahbration and the validation of 

b lied and accurately measured in Local G roup (d < 1 Mpc) and the PL relation should be rooted on distances measured with a 

■ ■ ■ ' in Local Volume (d < 10 Mpc) g alaxies (iFreedman et al.ll200Tl geometrical method such as the trigonometric parallaxes. This 

iTammann et al.]l200l iBono et alJiSl) . They obey to a Period- approach was recently adopted by Benedict et al. (2007) who 

Luminosity (PL) relation and are the most popular primary dis- provided parallaxes with a mean accuracy of 8% for a sample of 

tance indicators ('peas? 199?; 'Macri et al." '2006'; 'Fouque et al." "me Galactic Cepheids using the Fine Guidance Sensor (FGS) 

2007; Groeneweaen 2007; di Benedetto 2008; Groenewegen on boai'd of the Hubble Space Telescope (HST). A new revi- 

2008; Kervella et al. 2008; Kanbur et al. 200^ IMarengo et al.i sion of HIPPARCOS parallaxes for Galactic Cepheids (244 ob- 

2009: Scowcroftet al..,2009i) . In spite of these o utstanding ob- jects) has been rece ntly pr ovided by van Leeuwen ( 2007) and by 



servational and theoretical efforts (lMarconii2009l and references [van Leeuwen et^ The accuracy of the new measure- 

therein) the universality of optical and near-infrared (NIR) PL ments is on average a factor of two better than the old ones. 

However, the Cepheids with the most accurate HIPPARCOS par- 
allaxes are Polaris {a UMa) and the prototype S Cep. The accu- 



* Based on observations made with MPG/ESO 2.2m telescope at racy of the former one is 1.6% (van Leeuwen et al. 2002), while 
La Silla Observatory under proposal IDs: 75.D-0676, 60.A-9120 and the latter is similar to the accuracy of the FGS @HST parallax, 
multi-epoch, multi-band NIR data at SAAO. namely 5.2% versus 4.1% (iMerand et al.ll2005h . 
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The Baade-W esselink method (BW, iBaadd (Il926h : 
IWesselin^ (Il946h ) provides an independent empirical ap- 
proach to measure Cepheid absolute distances and it can be 
applied to variable stars. This method relies on two observables: 
the radial velocity (vv) and the variation of the angular diameter 
(ff). The latter parameter was historically substituted by the 
variation in color along the pulsation cycle. However, direct 
measurements of the Cep heid angular diam eter have been 
recently provided by Ker vella et al.l ( l2004bl) using the Very 
Large Telescope Interferometer (VLTI). In particular, the use 
of VLT INterferometer Commissioning Instrument (VINCI) 
gave the unique opportunity to provide angular diameter 
measurements along the pulsation cycle for seven Cepheids. 
These pioneer measurements encouraged a detailed comparison 
between theory and observations concerning the limb darkening 
and the atmosphere of variable stars (Marengo et al. 2003; 
iNardetto et al.l |2007|) . but the number of Cepheids for which 
these measurements are available is still very limited. 

To overcome the difficulties in the interferometric measure- 
ment of the angular diameter several variants of the BW method 
were suggested in the literature. Among them the methods 
based on the Surface-Brightness (SB) relations link variations 
in color with variations in angular diameters. This method can 
be applied to Cepheids for which accurate radial veloci t ies an d 
multi-wavelength l ight curv es are availab le. IStormetalJ(l2004 . 
iKervella et al]l20b4a) and .Groenewegenl(l2007h (hereafter G07) 
derived such relations, using an optical-NIR (V-K) color (IRSB), 
which gives the highest precision in the derived quantities. 

However, the most relevant limit of the currently adopted 
BW methods is the value of the projection factor (p-factor). 
This parameter links radial velocity changes to radius changes 
and still lacks firm theoretical and empirical constraints 
(INardetto et al.l l2009l) . It can be empirically estimated using 
Cepheids for which accurate interferometric angular diameter 
measurements, radial velocity curves and trigonometric paral- 
laxes are available. This approach was applied to S Cep by 
iMerand et al.l (l2005l) using the new optical interferometric mea- 
surements obtained with the CHARA Array, t he FGS@HST 
trigonometric parallax by ' Benedict et al.l (l2007 h and the radial 
velocities available in the literature. The /7-factor they found - 
p=1.27 + 0.06- agree s quite well with theoretical predictions by 
INardetto etaP (l2004 . More recently, G07 found that the use of 
a constant p-factor (p-1.27 + 0.05) for six Galactic Cepheids, 
with interferometrically measured angular diameter variations 
and known distances, agrees quite well with HST parallaxes. 
Moreover, he found that a s trong period depen dence of the 

p-factor (p 0.15 ■ logf , iGieren et all (l2005h ). could also 

be ruled out. However, a moderate period dependence (p ~ 
-0.03 • logP) as suggested either by Gieren et alj (I1993L 1 19971 
ll998l):lBames"et al (20C)3) ; Storm et al (2004), or more recently 

bv INardetto et all (l2009T ( p 0.08 • logP) is still consistent 

with currently available data (G07). 

In this investigation we plan to apply the BW method using 
the most recent calibration of the IRSB relation (Groenewegen 
2010, hereinafter GIO, in preparation), to estimate the distance 
of four metal-rich Galactic Cepheids. In particular, we plan 
to use new radial velocity measurements, new accurate NIR 
(J,H,K) l ight curves and V-band light curves available in the 
literature (lBerdnikovll992l) . Moreover, we discuss in §4 the iron 
abundance of the four target Cepheids using Fe i and Fe ii lines. 
§5 deals with the pulsation amplitude of metal-rich Cepheids, 
while in §6 we summarize current findings and outline future 
developments of this project. 
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Fig. 1. From left to right J,H, K-hand light curves for the four 
metal-rich selected Cepheids. In each panel are also plotted the 
intrinsic scatter (cr) of the fit with a cubic spline (red line) and 
the luminosity amplitude. 



2. Photometric and spectroscopic data 

The selected stars are the most m etal-rich Cepheids i n the larg e 
sample (113 objects) collected bv lAndrievskv et al] (l2002allbUcl 
2004). Positional and physical parameters for the selected 
Cepheids are listed in Table[T] 



2.1. Optical and NIR data 

In order to provide accurate estimates of their distances we 
plan to use the mo s t recent calibrat i on of t he IRSB relation 
(Bai-nes et al.' '1976'; 'L anev & Stobi^ 119951; iFouque & GierenI 
11997; Storm_et al. 2004 1 GIO). To accomplish this goal one of 
us (CDL) collected accurate multi-epoch, multi-band NIR data 
at SAAO for the four targets. The typical uncertainty of individ- 
ual phase points ranges from 0.005 to 0.007 mag for K < 6 mag, 
deteriorating to about 0.012 at K = 8.6 mag. Figure [T] shows the 
J, H and K-hand light curves for the target Cepheids. The in- 
trinsic scatter (cr) of the fit with a cubic spline (red line) and the 
luminosity amplitude are also labeled. The mean NIR {J,H,K) 
magnitudes have been estimated as a time average on magni- 
tudes and are listed in Table [T] The typical uncertainties are of 
the order of a few hundredths of a magnitude. 

For three out of the four targets (AV Sgr, UZ Set, 
VY Sgr) we have detailed optical light curves collected in 
the W 'ah-aven bands (Wa lraven et al. 1964; Lub & Pel 1975; Pel 
ll976l:lLub &. Pell 1977^ However, for the Baade-Wessehnk solu- 
tion we decided to use more recent V-band light curves available 
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Table 1. Intrinsic parameters, mean NIR magnitudes and mean radial velocities for the target Cepheids. 



NAME 


Q'(J2000)" 


SimOO)" 


logP 


ECB-V)" 


< J >' ±(t(J) 


<H >^ ±cr(H) 


< K >'■ ±cr{K) 


< 


< V,. >'' 




V340 Ara 


16 45 19 


-51 20 33 


1.32 


0.574 


7.382+0.011 


6.809+0.007 


6.619+0.008 


25-1-2 


-80.8+1.2 


59.1 


UZSct 


18 31 22 


-12 55 00 


1.17 


1.071 


7.502+0.049 


6.818+0.042 


6.564+0.045 


25-1-2 


40.2+0.5 


49.4 


AVSgr 


18 04 49 


-22 43 00 


1.19 


1.267 


6.909+0.040 


6.081+0.035 


5.758+0.033 


25-1-1 


19.3+1.8 


59.3 


VYSgr 


18 12 05 


-20 42 00 


1.13 


1.283 


7.174+0.069 


6.375+0.046 


6.068+0.040 


25-1-2 


16.0+1.8 


59.2 



" Cepheid coordinates: units of right ascension are hours, minutes, and seconds; units of declination are degrees, arcminutes and arcseconds. 
* Reddening according to FemieetaL Q995). 
^ Mean NIR magnitudes. 

Number of spectra collected with FEROS at the 2.2m MPG/ESO telescope. 

Mean radial velocity and velocity amplitude (km s"'). 



V340 ARA A= 59.1 <Vr>=-80.8 ±1.2 




° r UZ SCT A= 49.4 <Vr>= 40.2 ± 0.5 

V 20 - /^^"^ 
60 - ^^..V 



AVSGR A= 59.3 <Vr>= 1 9.3 ± 1 .8 




VYSGR A= 59.2 <Vr>= 1 6.0 ± 1 .8 

-20 - 
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Fig. 2. Radial velocity curves for the four selected Cepheids, 
measured by cross-correlating a line list (see text). The red aster- 
isks display the radial velocities based on the high S/N spectra 
adopted for measuring the iron abundance. The amplitude and 
the mean radial velocity (kms"') are labeled. 



in the literature (lBerdnikovll992h . since the accuracy improves if 
photometric and radial velocity data are collected close in time. 
The large sample of Cepheids with Walraven photometry will be 
used in Section 5, where we discuss the period-amplitude dia- 
gram. 



Table 2. New weak Fe ii line list added to the Fe i and Fe ii line 
list of ,Romani ello et al. (200l). 



J r A 1 


— 

Ion 


t^r 


— 

log 8.f 




Fen 


Z.5J 


-4.4j 




Fen 


TOO 

2.50 


1 Q C 




Fe II 


O 1 

z.ol 


-i.jo 


jlUU.OD 


Fe II 


O O 1 


-4. to 


^ 1 AT 


Fen 


O 1 

z.oi 


-i.yj 




Fen 


i.li 


-Z.Zi 


JZJ4.03 


Fe II 


D.IZ 


-L.LL 


jIjo.yA- 


Fe II 


Ton 


-4.25 




Fen 


i.ll 


^ 1 o 

-J. is 


3414. U/ 


Fen 


i.ll 


-J.j4 


j4zj.ZD 


Fe II 


i.l 


-i.ll 


jjJ4.;5j 


Fe II 


1 1/1 
J.z4 


-2. / J 


jyil.yjO 


Fen 


i.l 


/I O 1 

-4.51 


7135.02 


Fen 


6.21 


-2.60 


7301.56 


Fen 


3.89 


-3.68 


7310.22 


Fen 


3.89 


-3.36 


161231 


Fen 


3.15 


-5.19 


1111. n 


Fen 


3.20 


-5.24 


7801.24 


Fen 


5.91 


-2.94 


7841.39 


Fen 


3.90 


-3.72 


8250.34 


Fen 


5.22 


-3.29 


8264.72 


Fen 


6.81 


-2.08 


8324.96 


Fen 


6.22 


-2.75 


8330.59 


Fen 


6.22 


-2.35 


8490.08 


Fen 


9.74 


0.34 


8981.13 


Fen 


6.72 


-2.18 


8990.39 


Fen 


6.23 


-2.38 


9095.11 


Fen 


9.65 


0.27 


9122.92 


Fen 


9.85 


0.53 


9132.37 


Fen 


9.85 


0.51 


9187.16 


Fen 


9.70 


0.25 


9244.74 


Fen 


6.22 


-2.40 


9297.27 


Fen 


9.65 


0.41 


9464.88 


Fen 


6.09 


-2.54 


9572.60 


Fen 


5.82 


-2.87 


9843.19 


Fen 


6.14 


-2.60 


9849.74 


Fen 


6.73 


-2.30 


9956.31 


Fen 


5.48 


-2.95 



2.2. Spectroscopic data 

Together with photometric data we also secured for the tar- 
get Cepheids accurate multi-epoch, high-resolution (R~30,000) 

spectra covering the entire pulsational cycle with FEROS @2. 2m were collected between May and September 2005. The typi- 

MPG/ESO telescope. The goal was to provide individual radial cal exposure time is ~ 320s and the quality of the individual 

velocity measurements with an accuracy of 1 km/s for w two spectra is quite good (signal-to-noise ratio, S /N ~ 70) and al- 

dozen of epochs. The number of epochs is a strong requirement lowed us to reach the required accuracy in radial velocity mea- 

to reach with BW methods the nominal accuracy of 5% in dis- surements. Radial velocity measurements were performed using 

tance. The spectroscopic sample, indeed, includes more than 25 fitline a semi-inter active routine developed by one of us (PF, 

spectra per Cepheid (see column 9 in Table [T]i. These spectra lFran9ois et al.l (l2007h ). The algorithm adopted in fitline is 
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Fig. 3. High S/N spectra for the four selected Cepheids in a lim- 
ited wavelength range (5930-6000 A). Vertical dashed and dot- 
ted lines mark Fe i and Fe ii lines. In each panel is also plotted 
the S/N ratio of either individual (AV Sgr) or co-added spectra. 



based on a cross-correlation between the lines of each spectrum 
and a given line list. The reader int erested in a more deta i led de - 
scription of fitline is referred to lLemasle et al.1 (I2007ll2008h . 
Note that to pr operly deal with metal-ric h Cepheids the iron line 
list provided bv lRomaniello et al] (|2008[) was supplemented with 
three dozen of weak Fe ii lines (see Table|2]l. Figure |2] shows the 
individual radial velocity measurements and the fit with a cubic 
spline (red lines). The mean radial velocities (< v,. >) and the 
velocity amplitudes (Av,) based on the spline fit are also labeled 
(see also columns 10 and 11 of Table[T]i. 

The target Cepheids are among the most metal-rich Galactic 
Cepheids and, as expected, they are located close to the inner 
edge of the thin disk ( Pedicell i et al ]|200i. In order to pro- 
vide accurate iron abundance measurements, we also collected 
for each target 1-2 high S/N ratio spectra with FEROS@2.2m 
MPG/ESO telescope (see red asterisks in Figure |2] and Table 
[TJ. These spectra were collected by one of us (PF) at the end 
of March 2007. The typical exposure time is ~ 1800s and 
their S/N ratio ranges from 150 to 300. The abundance analy- 
sis of the Cepheids with two high S/N spectra was performed 
on the co-added spectrum, since they were collected close in 
time. Figure [3] shows the co-added spectra for the four selected 
Cepheids in the wavelength range between 5930 and 6000 A. 
The vertical dashed and dotted lines mark Fei (/l=5930.17, 
5934.66, 5956.70, 5976.78, 5983.69, 5984.79, 5987.05 A) and 
Fell (/1=5932.06, 5991.37 A) Hnes, respectively. The S/N ratio 
of either individual (AV Sgr) or co-added spectra is also labeled. 



3. BW method 

The BW method adopted to estimate distances and radii has 
already been discussed in previous papers by Groenewegen 
(Groenewegen 2004, 2007, 2008). In the following we briefly 
mention the key points of the method and the differences with 
the original approach. According to the definition of quasi- 
monochromatic surface brightness, we can write My - Sy + 
SxlogiR/Ro) - const. The absolute magnitude, and in turn the 
distance, can be found by differentiating this equation with re- 
spect to the pulsation phase, by multiplying the result for a color 
index -(V - K)()- and then by integrating over the entire pul- 
sation cycle. The radial velocity is tightly connected with the 
pulsation velocity and the curve can be integrated to obtain the 
radius variation as a function of time (phase); 



AR(t, 58) ■■ 



•Jta 



(V,{t) - Vy)dt 



(1) 



where P is the pulsation period, p the projection factor, 
Vy the systematic radial velocity and 66 accounts for a phase 
shift between the radial velocity curve and the angular diam- 
eter changes measured either interferometrically or via the SB 
relation. Eventually, we end up with the equation ©(?)= const. x 
[Rq + AR{t, 6ff)]/d where 0(f) is the angular diameter in mas, Rq 
is stellar radius in solar radii and d the distance in parsec. We 
measure the apparent radial velocity v, , i.e. the Doppler shift of 
absorption lines in the stellar atmosphere, projected along the 
line of sight and integrated over the stellar disk. To obtain the 

Silsational velocit y we use the p-fac tor (see column 2 in Table 
coming from .Nardetto et al.1 (l2009l) : 



p = 1.31 + 0.06 - 0.08 ± 0.05 ■ log(P) 



(2) 



Integration of the pulsational velocity over the entire pulsa- 
tion cycle provides an estimate of the linear radius variation. The 
angular diameter variation was derived using the IRSB method 
and the (V - K) color. The calibration adopted in this investiga- 
tion was derived by GIO: 



log 00 = 0.2692 (V - K)q + 0.5298 



(3) 



The above relation together with the Eq.[T] gives, once inte- 
grated over the pulsation cycle, individual Cepheid distances. 

3.1. Cepheid distances 

Figure |4] shows from left to right the application of the BW 
method to V340 Ara, UZ Set, AV Sgr and VY Sgr. For each 
Cepheid we plot the variation of the angular diameter against 
phase (bottom panels) and the change in angular diameter de- 
rived from the IRSB relation (see Eq. 13) against the change in 
radius (top panels) obtained by integration of the RV curve (see 
Eq. [TJ. The absolute magnitudes have been derived in a self- 
consistent way from the Fourier fit to the data. Note that the 
fit in the phase interval 0.8-1.0 is often poor. This limitation of 
the BW method was known before ( Storm et al.i,2004 ). It was 
argued that non-LTE eff'ects and an increase in the micro turbu- 
lence during these phases may change the atmospheric structure, 
and in turn h amper the use of a simple surface-brightness rela- 
tion (IBersier et al.lll997 ). Table [3] lists the estimated distances 
and radii and the comparison with previous estimates. The errors 
listed in columns 3 and 4 of Table [3] give the uncertainties (one 
cr) on distance and radius for the target Cepheids. They account 
not only for the errors in the fit, but also for the errors estimated 
running several Monte Carlo simulations by artificially chang- 
ing individual errors on V, K and RV measurements. In order 
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Fig. 4. From left to right: V340 Ara, UZ Set, AV Sgr and VY Sgr For eaeh star, the top panels show the linear-bisector fit to the 
angular diameter as a function of radial displacement. The bottom panels show the angular diameter against phase. Crosses mark 
data-points not included in the fit. 



to compare the above distances with similar estimates available 
in the literature we adopted the 7, H and /T-band PL relations 
for Galactic Cepheids recently provided by Fouque et al. ( 2007j). 
The mean of the three individual distances together with their 
errors for the target Cepheids are listed in column 5 of Table |3] 
The errors account for the uncertainty on the mean magnitudes 
and on the PL relations. The comparison indicates that distances 
based on the BW method and on NIR PL relation agree within 
Icr. Current uncertainties on the BW distances, by taking into 
account the errors in the fit and the errors estimated using the 
Monte Carlo simulations, range from 6% (UZ Set, AV Sgr, VY 
Sgr) to 14% (V340 Ara), while the uncertainties on the distances 
based on the empirical NIR PL relations range from 22% to 34%. 

The NIR PL relations provided by iFouque et all (l2007h do 
not account for a possible dependence of the NIR PL relation 
on the metal content. We are dealing with Cepheids character- 
ized by super-solar iron abundances. Therefore, we estimated 
the Cepheid distances using theoretical NIR PL relations in- 
cluding both metal-intermediate and metal-rich Cepheid mod- 
els (F. Caputo, private communication). Data listed in column 6 
of Table [3] show that distances based on these theoretical NIR 
PL relations also agree within Icr with the BW estimates. Note 
that the distances and their uncertainties based on theoretical PL 
relations are larger than the empirical ones, since the former re- 
lations were derived including all the fundamental models with 
abundances ranging from Z = 0.001 to Z = 0.04 and different 
helium contents (Bono et al. 2010, in preparation). 

The mean Cepheid radii listed in column 4 of Table [3] agree 
quite well with mean radii based on empirical Period-Radius 
(PR) relation for Galactic Cepheids recently provided by G07 
(see column 7 in Table O. Note that this PR relation i s based 
on six Galactic Cepheids with known distances (Benedi ct et alj 
l2007h and wi th interferometrically measured angular d iame- 
ter variations dKervellaet al.ll2004bt iMerand et al1l2005l) . The 
same outcome applies in the comparison with radii based on 
the predicted PR relation for Galactic Cepheids provided by 
[Petroni et al. (2003), assuming a solar chemical composition. 
Current radii agree on average within 10% (two cr) with pre- 
dicted and empirical PR relations. 



4. Cepheid iron abundances 

To measure the iron abundances of the t arget Cepheids we fol - 
lowed the same approach suggested by iLemasle et al.l (l2007h . 
This approach relies on three different steps: 

- Line list: The line list adopted bv iRomaniello et afl ( l2008l) . 
includes 275 Fe i lines and 37 Fe ii lines covering the FEROS 
spectral range. However, the target Cepheids are metal-rich 
and a good fraction of Fe ii are saturated. In order to provide 
robust estimates of intrinsic parameters and in turn accurate 
measurements of iron abundan ces, we supplemented the line 
list bv Roma niello et al.l ( l2008l) with 39 new weak Fe n lines 
(see Table|2]i. 

- Equivalent width: The measurement of the equivalent 
widths (EW) of the iron lines was performed using fit line. 
This code uses a Gaussian fit, which is defined by four pa- 
rameters: central wavelength, width, depth and continuum 
value of the individual lines. The initial value of the Gaussian 
parameters is fixed by randomly selecting the four parame- 
ters. Then, the "genetic" algorithm computes the;^'^ between 
the observed line and the expected Gaussian profile and com- 
putes the new set of Gaussian parameters among the 20 best 
fit solutions of the previous "generation" by applying ran- 
dom modifications of the values of the parameters ("muta- 
tion"). The algorithm after 100-200 "generations", gives the 
best fit Gaussian parameters (lowest x^) for each observed 
line. For the measurement of the iron abundance, we have 
selected only lines with equivalent widths between 10 and 
200 mA. The lower limit was chosen to be a safe compro- 
mise between the spectral characteristics and the need for 
weak lines for an optimal abundance determination. The up- 
per limit was fixed to avoid the saturated portion of the curve 
of growth. 

- Stellar Parameters: The determination of an accurate effec- 
tive temperature is a critical point in the abundance deter- 
mination. This requirement becomes even more important if 
we are dealing with variable stars, since the temperature es- 
timate has to refer to the pulsation phase at which the spec- 
trum was collected. In this analysis, the effective temperature 
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Table 3. Distances and radii according to the BW-analysis of the four selected Cepheids. 



Name 


P 


rf(pc) 




rf(pc)°, 


rf(pc)«P, 






V340 Ara 


1.205 


3890+126±547 


100.2±3.2+17.3 


3154+686 


3754±1156 


109.5+1.1 


115+1 


UZSct 


1.217 


3176+116±152 


82.1±3.0±3.5 


3117±871 


3540±1177 


86.4+1.1 


91+1 


AVSgr 


1.215 


2302± 72±123 


93.4±2.9±4.2 


3036±792 


3326±1063 


89.2+1.1 


94+1 


VY Sgr 


1.219 


2546± 53±134 


86.5±1.8+3.9 


2300±663 


2533± 849 


81.1 + 1.1 


86+1 



" Mean ///^'-distances based on empirical NIR PL relations (Fouqu e et al. 2007) . 

* Mean y/f/f-distances based on theoretical NIR PL relations (Caputo, private communication). 

Radius estimate according to the empirical PR relation by G07. 

Radius estimate according to the theoretical PR relation by^etroni^taL (2003). 



(T^eff) was estimated spectroscopically using the line depth 
ratios (LDR) method described in Kovtvukh & Gorlov^ 
(2000). This technique has the advantage of being indepen- 
dent of interstellar reddening and minimally dependent on 
metallicity. These uncertainties plague other methods like 
the integrated flux method or the color-temperature relations 
{Gray 1994; Krockenberger et al. 1998). This approach be- 
comes even more relevant for the target Cepheids, since they 
are located close to the edge of the inner Galactic disk, and 
therefore they are characterized by high reddening values 
(see column 5 in Table [T]). The estimated effective temper- 
atures together with their errors are listed in columns 2 of 
Table |4] The surface gravity (log g) and the microturbulent 
velocity (vt) were constrained by minimizing the log([Fe/H]) 
vs. EW slope (using the Fei abundance) and by impos- 
ing the ionization balance between Fe i and Fe ii (see col- 
umn 3 and 4 in Table |4|. These two procedures are tightly 
connected and require an iterative process. The initial val- 
ues for the microturbulent velocity and the surface grav- 
ity, were fixed using t ypical Cep heid values (vt=3 km 5"', 
log g=2. lAndrievskv et alj (l2002b)). For the ionization bal- 
ance, we assume that it was fulfilled when the difference be- 
tween [Fei/H] and [Feii/H] was smaller than the standard 
deviation on [Fen/H] (typically, cr[FeiH] ~ 0.08 - 0.1 dex). 
If this condition was satisfied by more than one value of 
log^, we checked which value satisfies also the ionization 
balance within the standard deviation on [Fe i/H] (typically, 
C[Fe/H] ~ 0.02 dex). The effective temperatures listed in col- 
umn 5 (reff(Model)) are the final best fit values of the at- 
mosphere models adopted in the iterative process. To deter- 
mine the errors on the microturbulent velocity and the sur- 
face gravity, we ran several iterations for each star, slightly 
modifying the values of these two intrinsic parameters that 
fulfill the requirements mentioned above. We have estimated 
that the intrinsic error on the microturbulent velocity is of the 
order of 0.1 km s ', while the intrinsic error on the surface 
gravity is of the order of 0.10 dex. 

Our final Fe i and Fe ii abundances, together with the adopted 
stellar parameters, are listed in Table |4] Data listed in column 8 
of this table indicate that the mean iron contents have weighted 
intrinsic uncertainties of the order of 0.20 dex, due to errors 
in the EW measurements, in the intrinsic parameters and in 
the number of unsaturated iron lines. This finding further sup- 
ports the difficulty in measuring the iron abundance of metal-rich 
Cepheids. The intrinsic accuracy can be certainly improved us- 
ing high S/N ratio, multi-epoch spectra of the same targets, since 
empirical evidence suggests that Cepheid e lemental abundances 
minimally depend on the p ulsation phase dLuck & Andrievskvl 
l2004HKovtyukh et alJl2005h . Finally, it is worth mentioning that 
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259 Goloctic Cepheids 



<[Fe/H]> = -Q.03 




Fig. 5. Top - Metallicity distribution of Galactic Cepheids. 
Bottom - Same as the top, but the metallicity distribution was 
smoothed using a Gaussian kernel with standard deviation equal 
to the metallicity uncertainty of individual Cepheids. The red 
line shows the Gaussian fit of the metallicity distribution. The 
mean and the cr are also labeled. 



current iron abund ances agree, within one cr, with the abun- 
dances provided bv lAndrievskv et aP (l2002bh using a similar ap- 
proach (see last column in Table|4]l. 



5. Period-Amplitude diagram 

The metal-rich regime of classical Cepheids has been only 
marginally investigated, since till a few years ago iron measure- 
ments for these objects were not available. As we have already 
mentioned in §2, Walraven VBLUW photometry is available for 
more than 160 Galactic Cepheids. The pulsation properties of 
these objects will be discussed in a forthcoming paper (Pedicelli 
et al. 2010). Interestingly, three out of the four target Cepheids 
belong to this sample, namely AV Sgr, UZ Set and VY Sgr For 
these Cepheids were secured at least 30 phase points in five- 
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Name 


reft(LDR) 


log 8 




r,ff(Model) 


[Fei/H] 


[Feii/H] 


<[Fe/H]> 


[Fe/H]^„,, 


V340 Ara 


5425+146 


0.2 


3.5 


5475 


0.38±0.17 


0.42±0.26 


+0.40+0.21 


+0.31+0.10 


UZSct 


4790+95 


0.6 


3.7 


4850 


0.35±0.17 


0.35±0.27 


+0.35+0.21 


+0.33+0.10 


AVSgr 


5407±53 


0.9 


4.2 


5450 


0.26±0.12 


0.28±0.12 


+0.27+0.12 


+0.34+0.10 


VY Sgr 


5268±146 


1.0 


3.2 


5400 


0.38±0.17 


0.36±0.26 


+0.35+0.21 


+0.26+0.10 



VBLUW bands that properly cover the entire pulsation cycle. 
The intrinsic accuracy of individual measurements is of the or- 
der of a few millimag. The uncertainty on the mean magnitudes, 
estimated using a fit with a cubic spline, is at most of the order of 
a few hundredths of a magnitude. In order to investigate the pos- 
sible correlation between luminosity amplitude and iron abun- 
dance, the Walraven BiyiViv-band photometric data were trans- 
formed into the Johnson B y,yy-band data using the relations: 



Vj = 6.886""'« - 2.5 Vw - 0.1916 (V - Bh 



(4) 



(B-V)j = 2.7947 (V-S)„,-1.2052 (V-By„+0M22 (V-fi)i^-0.0100(5) 

Accurate iron abundance measurements for Walraven 
Cepheids are available for 77 objects, while for 67 objects metal- 
licity estimat es are av ailable based on the Walraven metallicity 
index ( Pedicelli et al and references therein). This sam- 

ple was supplemented with 115 Cepheids for which accurate 
iron abundances exist in the literature, based on high resolution 
spectra (^nd rievsky et al. 2002a b.c, 2004; Lemasle et al. 2007; 
ISziladi et al.l2()07l:ILemasle et al. 2008; R omaniello et al.l2008l) . 
We ended up with a sample of 259 Galactic Cepheids and the top 
panel of Fig.|5]shows the metallicity distribution. In order to pro- 
vide accurate estimates of the mean metallicity and of the spread 
in metallicity of the Galactic disk, we ran a Gaussian kernel with 
a o" equal to the metallicity uncertainty of individual Cepheids. 
The metallicity distribution we obtained is plotted in the bot- 
tom panel of Fig. |5] and the red line shows the Gaussian fit. We 
found a solar mean metallicity (< [Fe/H] >~ -0.03) and a sigma 
of 0.18 dex. These values agree quite well wit h similar metal- 
licity distributions availab le in the literature dChiappini et al.l 
120011; ICescutti et"ai]|2007h . Our metallicity distribution is also 
asymmetric, and indeed the metal-poor tail is shallower than 
the metal-rich one. This finding agrees with predictions based 
on Galactic ch emical evolution models (H olmber g et al. 2007.; 
lYin et al.ll2009l) . 

To constrain the dependence of the pulsation properties on 
the metal content, we adopted the Bailey diagram, i.e. luminos- 
ity amplitude vs. pulsation period. Data plotted in Figures |6]and 
|7]show the V, B luminosity amplitudes of the Cepheid sample 
adopted by Pedicelli et al. (2009). The amplitudes are based 
either on Walraven photometry -transformed into the Johnson 
system (see §2.1)- or on data i n the Galactic Cepheid catalog 
provided by iFernie et"an (Il995h . To constrain the possible de- 
pendence of the pulsation amplitude on the metal content we 
selected three diff'erent sub-samples representative of the metal- 
poor ([Fe/H]<-0.30 dex) tail, of the metal-rich ([Fe/H]>0.13 
dex) tail and of the peak metallicity (-0.04 <[Fe/H]<-0.02 dex). 
The three samples roughly include two dozen of Cepheids. Data 
plotted in the top panels of Figures |6] and |7] display that both V 
and B-band amplitudes are not correlated with the metal con- 
tent. This feature does not agree with a well established the- 
oretical and empirical evidence of RR Lyrae stars (Bono et al. 
I2OO7I) . i.e. the prototype of low-mass, helium burning radial 
variables located inside the so-called Cepheid instability strip. 



Such a difference can be partially explained with the empir- 
ical circumstance that Galactic RR Lyrae variables cover al- 
most three dex in metal content, when moving from the halo 
to the bulge, while Galactic Cepheids in the disk only cover 
one dex. Moreover, the Hertzsprung progression causes a sys- 
tematic de crea se in the pulsati on amplitude s for periods across 
ten days (iPav ne-Gaposchkin 1951L 1 19541). Current empirical 
(Andreasen & Petersen 19871; IWelch e t al."1997:'Beaulieu il998t 
MoskaHketal. 2000) and theoretical (Bonoetal. 2000b|) evi- 
dence indicates that a decrease in metal content causes a sys- 
tematic drift in the center of the Hertzsprung progression to- 
ward shorter periods. This means a reshuffle in the pattern of 
the luminosity amplitudes as a function of the pulsation period 
when dealing with Cepheids with different metal abundances. 
However, the spread in the amplitudes is also present at peri- 
ods shorter and longer than the Hertzsprung progression. This 
indicates that the lack of a well defined correlation with metal- 
licity might also be due to the dependence of the topology of the 
instability strip on the ch emical composition (Pel & Lubl ll978l: 
iBono et al]|1999l l2000al) to evolutionary effects and to binarity 
(see Pedicelli et al. 2010, in preparation). 

As a final test concerning the metallicity distribution of 
Galactic Cepheids we split the sample into short (log P < 1.0) 
and long (log P > 1 .0) period objects. Data plotted in the bottom 
panels of Figures |6] and |7] show that metal-rich ([Fe/H]>-0.13 
dex) Cepheids are more frequent among long- than among short- 
period Cepheids (15% vs 5%) Cepheids. On the other hand, the 
metal-poor ([Fe/H]<-0.30 dex) Cepheids are equally distributed 
between the two groups. The evidence that the period distribu- 
tion of classical Cephe ids depends on the metal-content dates 
back to'Gascoigne' (1974, and references therein), who realized 
that the peak shifts toward shorter periods when moving from 
Galactic to Small Magellanic Cloud Cepheids. Subsequent evo- 
lutionary tracks for intermediate-mass stars showed that more 
metal-poor structures are characterized, at fixed mass, by blue 
loops that cover a larger temperature range. This means that the 
minimum Cepheid mass crossing the instability strip is smaller 
in metal-poor that in metal-rich systems. As a consequence a 
systematic drift in the period distribution of metal-rich Cepheids 
is expected. However, the period cut we adopted is well beyond 
the short period cutoff of Galactic Cepheids, thus suggesting that 
the paucity of metal-rich Cepheids among short-period Cepheids 
is probably due to an observational bias. 



6. Summary and conclusions 

We provided accurate BW distances and radii for four metal- 
rich Cepheids, namely V340 Ara, UZ Set, AV Sgr and VY Sgr. 
Current distance estimates, taken at face value, ag ree quite well 
with similar estimates based either on empirical (iFouque et alJ 
2007) or on theoretical NIR PL relations. However, the uncer- 
tainties affecting the BW distances, by summing in quadrature 
the errors in the fit and the errors estimated with the Monte Carlo 
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Fig. 6. Top - V-band amplitude versus period for a sample of 
259 Galactic Cepheids. Triangles and circles display spectro- 
scopic and photometric metallicities, respectively. Red and green 
dots mark metal-rich ([Fe/H]> 0.13) and metal-poor ([Fe/H]< - 
0.30) Cepheids, while the blue ones show Cepheids with iron 
abundances ranging from -0.04 to -0.02 dex. The selected metal- 
rich Cepheids are marked with a cross, while the pluses mark 
first overtone pulsators. Bottom - same as the top, but the ver- 
tical dashed line splits short- (log(P) < 1.0) and long-period 
(log(P) > 1.0) Cepheids. The fractions of metal-poor and metal- 
rich Cepheids are also labeled. 



simulations, are on average a factor of 3-4 smaller than for dis- 
tances based on predicted and empirical NIR PL relation. The 
same outcome applies to the mean Cepheid radii, but the uncer- 
tainties on the BW radii, by summing in quadrature errors on 
the fit and errors from Monte Carlo simulations, are on average 
a factor of two larger than for radii based either on the empir- 
ical or on the theo retical PR relation provided by GOV and by 
iPetroni et al.l (12003 ), respectively. 

We also collected high-resolution, high signal-to-noise ratio 
spectra to measure the iron abundances of the target Cepheids. 
Special attention was paid to provide accurate estimates of in- 
trinsic parameters (effective temperature, surface gravity, mi- 
croturbulent velocity) directly from observed spectra. We per- 
formed detailed measurements of iron abundances using large 
samples of Fe i and Fe ii lines. Current abundances indicate that 
selected Cepheids are super met al-rich and agree, within Icr, 
with iron abundances provided bv lAndrievskv et al] (l2002bl) us- 
ing a similar approach. 

We adopted a sample of 259 Galactic Cepheids for which 
are available either spectroscopic iron measurements or metal- 
licity estimates based on the Wakaven metallicity index 
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Fig. 7. Same as Figure|6l but for the B-band amplitudes. 



(iPedicelli et al.ll2008l) . We found that classical Cepheids do not 
seem to show in the Bailey diagram (luminosity amplitude vs 
pulsation period), in contrast with low-mass helium burning RR 
Lyrae stars, a clear correlation between luminosity amplitude 
and metallicity. The lack of such a correlation might be the con- 
sequence of the Hertzsprung progression. We also found that a 
good fraction of metal-rich ([Fe/H]>0.13 dex) Cepheids are lo- 
cated among long-period (logP > 1.0) variables. However, for 
the moment this can only be considered as circumstantial evi- 
dence, since the current sample is probably affected by selec- 
tion bias. Metal-rich Cepheids are located in the inner disk and 
are typically affected by large extinctions. The selected Cepheids 
have Galactocentric distances smaller than 6.5 kpc and their red- 
dening ranges from 0.6 to 1.3 mag. Detailed analysis concerning 
the pulsation properties of metal-rich Cepheids best awaits more 
complete samples. However, the game is worth the candle, since 
classical Cepheids are excellent tracers of young stellar popu- 
lations. Their pulsation properties and their radial distribution 
across the inner Galactic disk and the bar can provide robust 
constrain ts on the bar-driven formation scenario ( van Loon et al.l 
2003; De battista et ani2004t IZoccah et al.ll2006l) on short (10^ 
100 Myr) timescales. 
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